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Abstract 

The superfluid-Mott insulator transition of spin-2 boson atoms with repulsive interaction in an 
optical lattice in a magnetic field is presented. By using the mean field theory, Mott ground states 
and phase diagrams of superfluid-Mott insulator transition at zero temperature are revealed. Applied 
magnetic field leads to some phase boundaries splitting. For all the initial Mott ground states con- 
taining multiple spin components, different spin components take on different phase boundaries. It is 
found that in this system the phase boundaries with different magnetization can be moved in different 
ways by only changing the intensity of the applied magnetic field. 
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I. INTRODUCTION 



Recent remarkable experiments 



□ □ Q □ 



on the superfluid (SF) to Mott insulator (MI) 



transition in a system of ultracold atoms in an optical lattice open intriguing prospects for 
studying many-body phenomena, associated with strongly correlated systems in a highly con- 
trollable environment. The optical lattices — arrays of microscopic potentials induced by 
the ac stark effect of interfering laser beams — provide ideal conditions for the study of the 
laser cooling and the quantum phase transitions of the confined cold atoms. The dynamics 
of the confined atoms in optical lattices is adequately described by the Bose-Hubbard model 
[3, sj], which predicts SF-MI transition at low temperature with increasing the ratio of the 
on-site interaction to the hopping matrix element. Besides many experimental efforts made to 
realize SF-MI transition, a large number of theoretical studies have appeared |l, 0, 13 • In 
reference jsj an appropriate mean-field approximation was developed for the Hamiltonian of 
spinless or polarized bosons in an optical lattice, and in order to describe the zero-temperature 
phase transition from the superfluid to the Mott-insulating phase, the phase diagrams were 
calculated. 

Since optical traps [ill Q, Q liberate the spin degrees of freedom and make possible 
condensation of spinor bosons, extensive interests have been stimulated in the study of multi- 



component spinor BEC. Th e q uantum phase transition in spinor BEC, as well as a variety of 



other novel phenomena 




m nan 



13, 



, 3] were well studied. Subsequently, inspired by these 
works, the theoretical researches about the SF-MI transition of the spinor bosons trapped in 
an optical lattice arise. Most recently, Demler and Zhou have studied spin-1 Bose atoms 
in an optical lattice and obtained several unique properties. Tsuchiya et al. j^j], Hou and 
Ge [21] have investigated the spin-1 and spin-2 bosons in an optical lattice with the mean- 
field approximation method and obtained the phase diagrams showing a transition from Mott 
insulator to superfluid, respectively. 



On the ot 
interests 17, 



rer hand, the 



18, 



22 



23 



response to external magnetic field of BEC is also a topic with 



24, 



251 ]. The experimentalists have concentrated on investigating the 



systems in an applied magnetic field because the phase transition can be tuned by adjusting the 
magnetic field rather than changing the samples measured 2£ , 27. 0, 3 . Thus the theoretical 
study in this aspect is necessary. So far, Ueda and Koashi 



17 



18| have discussed the magnetic 



response of spin-1 and spin-2 BEC in a mesoscopic regime; as to the system in an optical 
lattice, Svicainsky and Chui have studied the spin-1 bosons in a magnetic field and shown 
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some effects induced by the magnetic field. However, the spin-2 case in a magnetic field has 
not been discussed yet. What about the ground states, quantum phase transition and the 
influence of the magnetic field on this system? This is our mission in the present letter. First, 
ignoring the hopping term of the Hamiltonian, we get the site-independent Hamiltonian, its 
energy eigenvalues and the Mott ground states for different cases. Then, applying mean-field 
approximation and regarding the hopping term as a perturbation, we perform the calculations 
in second-order and draw the phase diagrams. The response of the phase diagrams to the 
applied magnetic field is qualitatively analyzed subsequently. Finally, we give some remarks 
and the conclusion. 



II. THE MODEL 



We consider a dilute gas of boson atoms with hyperfine spin F = 2, such as 23 Na, 87 Rb or 85 Rb 



subject to an external magnetic field in an optical lattice. Based on Ref. [2l(, the Hamiltonian 
of spin-2 bosons with repulsive interaction in an optical lattice including a magnetic field term 
can be written in the second-quantized form: 

H = H A + H B , (1) 
Ha = J d ^^i • W a + V(r)^ a - fc^ a + |*^^ Q 

+f E(*i(^A) 2 +^i^^^ (2) 

i 

H B = -Hb9 J drtft ( B • F)^^, (3) 

where V(r) = Vo(sin 2 kx + sin 2 ky + sin 2 kz) with k the wave vector of the laser light and Vq a 
tunable amplitude characterizing an optical lattice, B is a uniform magnetic field independent 
of time, fiB is the Bohr magneton and g is the Lande factor of an atom, M is the atomic 
mass, \l/+2, ^-2 are the five- component field operators corresponding to the sublevels nip = 
+2, —2 of the hyperfine state F = 2, \x is the chemical potential, cq = 47rft, 2 (3a4 + 4a2)/7m, 
c\ = A-Kh 2 (di — 02) /7m and 62 = 47rfr 2 (3a4 — 10a2 + 7ao)/7m are parameters related to s-wave 
scattering lengths ao, ci2 and 0,4 of the two colliding bosons with total angular momenta 0, 2 
and 4, < 2a;2a'|00 > and < 00|2/?; 2/3' > are Clebsch-Gordan coefficients. F a (a = x,y,z) are 
5x5 spin matrices obeying usual momentum commutation relations [F a , Fp] = e a p 1 F 1 . 
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For simplicity, we assume a uniform magnetic field applied along the ^-direction, and it 
is weak enough to ignore the quadratic Zeeman effect. Expanding the field operators in the 
Wannier basis and keeping only the lowest vibrational states, = J2ibi a w(r — i-j), Eq. (JTJ) 
reduces to the generalized Bose-Hubbard Hamiltonian 

H = -t £ 6fe-/iE^+|E^(^- 1 ) + |E(F, 2 -6^ + ^E^--E^- ( 4 ) 

<i,j> i i i i i 

where Fj = bl a F a(3 b if3 , hi = £ Q b\ a b ia , t = -J drw*(r)(-h 2 V 2 /2m + V (r))w j(t) is the hopping 
matrix element between adjacent sites i and j, p = pj dr\wi(r)\ 2 + J drw*(r)(— h 2 V 2 /2m + 
V(v))wi(r) describes the effective chemical potential, and Cj = C{ j dr\wi(r)\ A is on-site inter- 
atom interaction, where the Hubbard approximation has been used to treat the multi-center 
integral as a single-center one. S i+ = Sj__ = (b\ ) 2 /2 — + ^2^1-2 creates a spin-singlet 

"pair" when applied to the vacuum, and we further assume p = gp-BB > in the following 
discussion. 

In the limit Co/t — ► oo, the hopping term can be neglected, so the Hamiltonian is reduced 
to a diagonal matrix with respect to sites. Then the single-site Hamiltonian is 

h = -ph + |n(n - I) + |(F 2 - Qh) + ^S + S^ - P F Z . (5) 

For the sake of studying the quantum transition, the mean-field approximation 9] is used 
and the hopping term is considered as a perturbation. Introducing the superfluid order pa- 
rameter (p a =< bi a >= y/n s f( a (n s f is the superfluid density and ( a is a normalized spinor 
= 1), we decouple the hopping term as b\ a bj a « {4>ab\ a + 4>abja) — 4>* a 4>a, then, the total 
hopping term becomes the product of a site-independent term and the total number of the sites. 
As a result, we only consider a single site because the Hamiltonian of every site is identical 
in the homogenous case, and the Hamiltonian is represented by a site-independent effective 
Hamiltonian multiplied by the total number of sites. The Hamiltonian of a single site reads 

h = h + hi, (6) 
hi = zt(cj) a bi + (p* a b a - (p* a (pa), (7) 

where hi is the mean-field version of the hopping Hamiltonian in the single site, and z is the 
number of the nearest-neighbor sites. When the ratio c /t is very large, hi is considered as a 
perturbation term. 
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III. THE ENERGY EIGENVALUES AND THE MOTT GROUND STATES 



Before perturbative calculations, we solve the equation 

h ip = e (0 V (8) 

to get the eigenvalues and eigenstates of h . In Eq. (JSJ), the operaters S + and S_ satisfy the 
SU(1, 1) commutation relations, namely [S z , S±] = ±S±, [S + , S-] = —2S Z together with S z = 

~ 2 » 

(2n + 5)/4, and as a consequence, the Casimir operator S = — S+S- + S 2 — S z commutes with 
S± and S z jl?! 3|- The eigenvalues of the mutual eigenstates for S and S z are {S{S — 1), S z } 
with S = (2n + 5)/4 (no = 0, 1,2, ...) and S z = S + n s (n s = 0, 1,2, ...), which guarantees 
that S + S- = S 2 — S z — S 2 is positive semidefinite. The new quantum numbers n s and n are 
introduced and the operator S+ raises n s by one while the relation n = 2n s + no holds, where n 
is the total number of bosons in a single site. The spin operator F and the magnetic quantum 
number operator F z commute with the operators S±, and then, the eigenstates ip are denoted 
as |n , n s , F, F z ; A > where A labels orthonormal degenerate states. The energy eigenvalue is 
given by 

e (°) = -/xn + jn{n - 1) + j[F(F + 1) - 6n] + -j-n,{n - n s + -) - pF z , (9) 

So Mott states can be expressed as Hi \ n o, n s, F, F z ; A >j, where i is lattice site index. For 
the homogenous case, the zeroth order total energy is E {0) = £;£ (0) = We®, where N t is the 
number of lattice sites. 

The energy eigenstates \n , n s , F, F z ; A > can be represented as [2 

(F^^r-p^o^ (io) 



where 



4 2)t = ^[( & o) 2 - ^\bU + 2bibi 2 i (ii) 



1 i2V2bibi-V3(bi) 2 }, (12) 



if t = _i=[V2(6 ) 3 - 3V2&I&0&U + 3v / 3(6 t 1 )V 2 + 3V3fo|(&-i) 2 - 6V26J&J&1.J, (13) 
^[(6 t 1 ) 3 -v / 646 t 1 6 + 2(6 2 )V 1 ]. (14) 



F(n s =o) is the projection onto the subspace with n s = 0; rii 2 , ^20, ^22, ^30 = 0, 1, 2, 00, n 33 = 
0, 1, and AF = 0, 1, 2F that are related to {n , n s , F, F 2 } by 

n = n 12 + 2n 22 + 3n 30 + 3n 33 , (15) 



n s = n 20 , (16) 
F = 2n 12 + 2n 22 + 3n 33 , (17) 
F z = F-AF. (18) 

From Eq. we see that the minimum energy states always satisfy F z = F when p > 0. 
Thus, the problem of finding the ground states reduces to minimizing the function: 

CO)/ TP \ , C / ,\ , C lr/i7 2p — Ci 2 (ci — 2p) 2 
e w (F,n s ) = -/m + -n(n-l) + -[(F ^—) - 6n ^ ] 

+ 2 -fn s {n-n s + \)- (19) 

The ground states of Hq depend on the relation among ci, c 2 ,p and \x. When n — 1, the ground 
state is |1, 0, 2, 2; A >; when n > 2, there are four sorts of ground states classified by different 
sign combinations of c\ and c 2 . Here the classification of the Mott ground states differs from 
that of Ref. 3], which is labelled by ferromagnetic, antiferromagnetic and cyclic phases. In 
general, for p > 0, when c\ < 0, c 2 > (ferromagnetic case), the third and the fourth term in 
Eq. (fTTJj) possess the minimal values simultaneously, but when c\ < 0, C2 < (ferromagnetic 
or antiferromagnetic case), c% > 0, c 2 > (ferromagnetic or cyclic case) and ci > 0, c 2 > 
(ferromagnetic or antiferromagnetic case), the situation becomes complicated, since sometimes 
there exists the competition between contributions of total spin and the singlet "pairs" to 
eigenenergy. We list all Mott ground states in detail in appendix. 

IV. PHASE DIAGRAMS OF SUPERFLUID-MOTT INSULATOR TRANSITION 

From Section III, we know that F = F z always holds in the ground states, while in the 
excited states this is not the case. As a result of the applied magnetic field, the state degeneracy 
from the different magnetic quantum numbers F z |2l| is lifted. More rich phase diagrams are 
expected than those in the absence of the magnetic field. 

To depict the phase diagrams, we consider the hopping term as the perturbative one and 
calculate the first- and second-order corrections to the ground energy, which are expressed as 

=< glhlg >= zt J2 0«, a = -2, 2, (20) 



a 



< g\hi\m > \ 2 s^sr^ z2 t 2 \ < dl^a + b^m > 



g — (o) (o) ~ 2^ 2^ (o) (o) , a— z, yzi) 



e 9 Zlv (0) (0) ^ (0) (0) 

n^g Eg Era n^g a Eg Em 
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Here \g > denotes the ground state discussed in the appendix, and {\m >} represent excited 
states expressed as a cluster of quantum numbers including n , n s , F, F z , A. We can calculate all 
the nonzero matrix elements of < g\b a + b' a \m >. Therefore, second-order perturbation theory 
gives the form of the modified ground energy as 

e g = ef + e« + ef = ef + zt£ 4*(n, £, c h c 2 ,P)<f>* a <f>*, a = -2, 2, (22) 

a 

where A a (n, jl, c , C\, c 2 ,p) is related to the first- and second-order corrections of the spin compo- 
nent with magnetic quantum number a to the zeroth-order ground energy. It depends on system 
parameters n, jl, cq, c±, c 2 and p, where jl = fi/zt,c.Q = c^jzt,C\ = c\jzt,c 2 = c 2 /zt,p = p/zt 
are dimensionless. Minimizing the ground energy function Eq. (|22)1. we find that (f) a = 
when A a (n, p,,c ,ci,C2,p) > and (fi a ^ when A a (n, p,, c , ci, c 2 ,p) < 0. This means that 
A a (n, jl, Co, ci, c 2 ,p) = signifies the boundary between the superfluid and the Mott insulator 
phases of the spin component with magnetic quantum number a. 

Using the perturbation theory, we can analytically determined the phase diagrams Fig. 
for different cases. The phase diagrams indicate that there exists a phase transition from Mott 
insulator with integer number bosons to superfluid when the ratio co/t is decreased to a critical 
value. In the zeroth-order, i.e., neglecting the hopping term, the ground state is Mott state in 
which the occupation number per site is pinned at integer n = 1, 2, corresponding to a com- 
mensurate filling of the lattice. Different ground states may contain different spin components. 
For example, there is only spin component with Zeeman level m = 2 when occupation number 
per site n = 1; spin components with m = 0, ±1,±2 for Mott state rL(^o |0 an d spin 

* (2) 

components with m — 0, 1, 2 for PL (^2 |0 >)»■ F° r the initial Mott ground state including only 
one spin component, one superfluid component occurs when lowering the ratio Co/t, such as the 
case n = 1, 2, 3 in Fig. [J n — 1, 2 in Fig. |21 and n — 1 in Fig. 010 for all the initial Mott ground 
states containing multiple spin components, when lowering the ratio co/t, multiple superfluid 
components appear, and the phase boundaries between superfluid and Mott insulator phase 
for different spin components are distinct, for instance, n = 3 in Fig. |21 n = 2, 3 in Fig. EE 
After analyzing the phase diagrams, we find that the position of phase boundary is related to 
average occupation number of spin component in the initial Mott ground state, i.e., the larger 
the average occupation number of spin component per site is, the easier the transition from 
Mott insulator to superfluid phase. We also find that some boundaries between superfluid and 
Mott insulator phases with multi-spin components, such as n = 2 in Fig. 13151 and n = 3 in Fig. 
13141 will turn to be identical when the magnetic field vanishes. We can draw the conclusion 
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that the applied magnetic field results in some phase boundaries splitting. 

Furthermore, the phase diagrams Fig. Efa) and Fig. Efb) for different intensities of the 
applied field are drawn. The position of some phase boundaries is related to the intensity of 
the applied magnetic field. In Fig. Efa), when the applied magnetic field increases, for the 
same Mott ground state containing only one spin component with Zeeman level m — 2, the 
phase diagrams will shift along the direction with chemical potential decreasing. In Fig. E^b), 
for the same Mott ground state containing spin components corresponding to Zeeman levels 
m — 0, ±1, ±2, when the magnetic field increases, the phase boundaries of spin components 
with Zeeman levels m = ±1, ±2 move, but the phase boundary of spin component with Zeeman 
level m = keeps invariant. Moreover, one can see that for positive and negative Zeeman levels, 
the phase boundaries between SF and MI will move in opposite directions, and for the spin 
component with positive Zeeman level, the transition from MI to SF becomes easier when the 
applied magnetic field increases. 

V. REMARKS AND CONCLUSION 

For simplicity, in this paper we assume the uniform magnetic field is applied along the z- 
direction, and p is positive. In fact, the analysis of the Mott ground states for p > is enough 
since the sign of p does not alter the physics. When p < 0, the Mott ground states satisfy 
F z = —F, and have the same form as those in the case of p > 0. We find that the influence of 
magnetic field on the phase diagrams is manifold, and the boundaries between SF and MI are 
essentially dependent on the magnetic properties of the ground states. This work only discuss 
the phase diagrams corresponding to different magnetic fields with the same Mott ground state, 
i.e., the continuous change of the phase boundaries. The case that Mott ground states change 
when the magnetic field increases, i.e., the sudden jump of the phase boundaries, is beyond the 
scope of the present paper and is not shown. For instance, when c\ > 0, c<i > 0, if p increases 
from to the value large enough to satisfy p/(2n + 1) > c 1; the magnetization F can jump 
from the minimum to the maximum one. They are the subject of future study. In addition, 
it is worth to note that phase diagrams of the zero magnetic field cannot be derived by taking 
p = simply, since our derivation is based on the degeneracy lifting. 

In conclusion, we have investigated the quantum phase transition from Mott insulator to 
superfluid phase of spin-2 cold bosons with repulsive interaction in optical lattices under the 
influence of a uniform magnetic field at zero temperature. The phase diagrams show that 
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the system undergoes a phase transition from Mott insulator with integer number bosons at 
each site to superfTuid phase when the ratio co/t is decreased to a critical value. Different Mott 
ground states may contain different spin components. The position of phase boundary is related 
to average occupation number of spin component in the initial Mott ground state. For the initial 
Mott ground state including only one spin component, one superfluid component appears when 
lowering the ratio Co/t. For some Mott ground states with multiple spin components, the applied 
magnetic field leads to the splitting of the phase boundaries, so that the phase boundaries 
between superfluid and Mott-insulator phase for different spin components are distinct in all 
ground states. In particular, we draw the phase diagrams corresponding to different intensities 
magnetic field for the initial Mott ground state containing one-spin and multi-spin components. 
They qualitatively show the way of the phase boundaries' moving with the intensity of the 
applied magnetic field. It is found that the phase boundaries can be moved by only changing 
the intensity of the applied magnetic field. For the spin component with positive Zeeman level, 
the larger intensity of the magnetic field is, the easier the transition from MI to SF happens. 
These theoretical results are expected to be practically helpful to the experimental study of the 
field-tuned SF-MI transition of bose atoms with hyperfine spin in an optical lattice. 
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APPENDIX: LIST OF THE MOTT GROUND STATES 

1. C\ < 0, C2 > 0. Because (2p — c\)/(2c\) < 0, the third term and the fourth term in Eq. 
(I19|) have the minimal values at the same time when F = 2n and n s = 0. Hence, the 
ground state is \n, 0, 2n, 2n; A >. 

2. ci < 0, c 2 < 0. The relation {2p — c\)/{2c\) < stands, and the competition between 
contributions of total spin and the singlet "pair" happens in the premise of ~ + 2n s = n. 
If we skip over the fact for the moment that the singlet "pair" number n s is an integer, 
then the condition to minimize the energy function is given by 



Because the singlet "pair" number must be an integer, we write n s in terms of the closest 
integer number n° s and the decimal part, i.e. n s = n° s + a, where the number a satisfies 
—1/2 < a < 1/2, which can be rewritten as, 

1 10c 1 (4n + l)-c 2 (2n + 3)-20c 1 p 1 
~ 2 < 80 Cl - 4c 2 < Us + 2' ( j 

(a) Because the singlet "pair" number is not negative, n s must be zero for n° < 0. 
Hence, the ground state is \n, 0, 2n, 2n; A >. 

(b) For < n° s < |, the eigenenergy is lower when the singlet "pair" number is n® than 
any other integer. So the ground state is |n — 2n°, 2n — 4n°, 2n — 4n°; A >. 

(c) For n° s > |, the singlet "pair" number takes the highest value as it can. (i)When n 
is even, the ground state is |0, n/2,0, 0; A >; (ii)when n is odd, the ground state is 
|l,(n-l)/2,2,2;A>. 

3. Ci > 0, c 2 > 0. For convenience, we introduce a new parameter 

F 1 = ^f^ = F + a, (A.3) 
2ci 

where F is the integer number which is the closest integer to Fi. In Eq. (|19|) . when 
n s = the fourth term has the minimal value, and if total spin F takes the appropriate 
integer the third term possesses the minimal one. 

(a) When F < 0, i.e., C\ > p, the energy eigenvalue has its minimum if F = and n s = 
0. However, we must notice some special cases because there have some forbidden 
values [l?], W^, that is, F = l,2,5,2n — 1 are not allowed when n = 3k(k G Z), 
and F = 0, 1, 3, 2hq — 1 are forbidden when no = 3k ± l(k G Z). 

i. For n = 3k(k G Z), |n, 0, 0, 0; A > with n s = and F = is the ground state. 

ii. For n = 3k — l(k G Z), the state with n s = and F = is not allowed 
simultaneously. If n s = 0, the lowest allowed value of total spin F is 2. On 
the other hand, F = is not forbidden when n is 3k(k G Z). n s = 1 is the 
lowest value satisfying the condition due to n = 2n s + n . So the state with 
n s = 1 and F = is a possible ground state. Comparing both eigenenergies 
for the two cases F = 0, n s = 1 and F = 2,n s = 0, we get the ground state, 
(i)|n - 2, 1,0,0; A > for c 2 < (15ci - 10p)/(2n+ 1) and (ii) |ri,0,2,2;A > for 
c 2 > (15c!-10p)/(2n+l). 
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iii. For n = 3k+l(k G Z), the competition exists between contributions of total spin 
and the singlet "pair" to eigenenergy. But if F = 0, n s is at least 2. Therefore, 
the ground state is (i)|n - 4, 2,0,0; A > for c 2 < (15ci - 10p)/2(2n - 1); (ii) 
|n, 0, 2, 2; A > for c 2 > (15c! - 10p)/2(2n - 1). 

(b) When F > 2n, i.e., c\ < p/(2n + 1), the ground state is |n, 0,2n, 2n; A > with 
F = 2n and n s = 0. 

(c) When < F < 2n, i.e., pj (2n + 1) < c\ < p, the eigenenergy is thus lower when F 
is closer to F\ and when n s is smaller. Except for the case of F taking the forbidden 
values of F, the ground state is \n, 0, F , F ; A >. When F equals to the forbidden 
values of F, F may take the allowed integer next-nearest to Fx, i.e., Fq ± 1 or Fq ±2, 
or at the cost of increasing n s to 1 or 2, since any of the three values 0, 1, 2 of n mod 
3 is realized by setting n s as 0,1, or 2 owing to the relation n = 2n s + no- Whether or 
not the states \n, 0, F ± 1, F ± 1; A >, |n, 0, F ± 2, F ± 2; A >, |n — 2, 1, F , F ; A >, 
and \n — 4, 2, F , F ; A > can be the lowest-energy state depends on the ratio c 2 /ci. 

4. c\ > 0, c 2 < 0. The eigenenergy is the lowest if F = F and n s at its highest value. 
However, these two choices are not always satisfied simultaneously. 

(a) Fq < 0, i.e., ci > p. 

i. When n is even, |0, n/2, 0, 0; A >is the ground state. 

ii. When n is odd, n s has the highest value (n — l)/2. But F is not zero when n s = 
(n — l)/2. Alternatively, there is another case that F = and ri s = (n — 3)/2. 
Hence, the ground state is (i) |l,(n — 1)/2,2,2;A > for ci < (7|c 2 | + 10p)/15; 
(ii) |3, (n - 3)/2, 0, 0; A > for c x > (7|c 2 | + 10p)/15. 

(b) Fq > 0, i.e., c\ < p. Note that n° s is shown in Eq. (jA.2jl . similar to the analysis of 
case 2., it is divided into three cases. 

i. For n° < 0, the ground state is \n, 0, 2n, 2n; A >. 

ii. For < n° < |, the ground state is |n — 2n°, n°, 2n — 4n°, 2n — 4n°; A >. 

iii. For n® > |, (i)when n is even, |0, n/2, 0, 0; A > is the ground state; (ii)when n 
is odd, 1 1, (n - l)/2, 2, 2; A > does. 
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FIG. 1: The phase diagram of Bose-Hubbard Hamiltonian obtained from second-order perturbation 
theory with solid lines for c\ = — Q.lzt, C2 = 0.1 zt and p = 0.2zt. The dashed lines indicate the 
zeroth-order phase diagram. 
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FIG. 2: The same as in Fig^ but for c\ = — 0.02zi, C2 = — 0.25zt and p = O.Olzt with solid lines. For 
n = 3, the interior line is the phase boundary of spin component with Zeeman level m = 2; the middle 
line the phase boundary of spin component with Zeeman level m = —2, and the external triple lines, 
which are too close to be distinguished, the phase boundaries of spin components with Zeeman levels 
m = 0, ±1. The inset shows an expansion of the region labelled by A. 
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FIG. 3: The same as in FigQ but for c\ = 0.22zi, C2 = O.lzt and p = 0.05zt with solid lines. 
When n = 2, the region labelled by B is enlarged in the left inset; the five lines represent the phase 
boundaries of spin components with Zeeman levels m = —2,-1,0,1,2 respectively (from interior to 
external). For n = 3, the five lines are too close to be distinguished, so the region labelled by C 
is enlarged in the right inset, in which the five lines the phase boundaries of spin components with 
Zeeman levels m = 0, ±1, ±2. 
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FIG. 4: The same as in Fig^ but for c\ = 0.18zt, c% = l.Ozt and p = 0.05zt with solid lines. For 
n = 2, the three lines represent the phase boundaries of spin components with Zeeman level m = 1,0,2 
respectively (from interior to external). When n = 3, the five lines are too close to be distinguished, 
so the region labelled by D is enlarged in the inset, in which the external line is the phase boundary 
of spin component with Zeeman level m = 2, and the other four lines the phase boundaries of spin 
components with Zeeman levels m = 0, ±1, —2. 
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FIG. 5: The same as in Fig^ but for c\ = 0.02zi, C2 = —Q.lzt and p = 0.008zt with solid lines. When 
n = 2, the five lines are too close to be distinguished, so the region labelled by E is extended in the left 
inset, in which the interior line is the phase boundary of spin component with Zeeman level m = — 1; 
the external line the phase boundary of spin component with Zeeman level m = 1; the three middle 
lines the phase boundaries of spin components with Zeeman levels m = 0, ±2. For n = 3, the exterior 
line is the phase boundary of spin component with Zeeman level m = 2; the middle line the phase 
boundary of spin component with Zeeman level m = —2, and the external triple lines, which are too 
close to be distinguished, the phase boundaries of spin components with Zeeman levels m = 0, ±1. 
The right inset shows an expansion of the region labelled by F. 
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FIG. 6: The phase diagrams of Bose-Hubbard Hamiltonian obtained from second-order perturbation 
theory for different p. In (a), c\ = — O.lzt, c% = O.lzt, with p = 0.2zt (solid lines) and p = O.Szt 
(dashed lines) respectively. In (b), when n = 2, c\ = 0.22zt, C2 = O.lzt, with p = 0.05zt (solid lines) 
and p = 0.082ii (dashed lines) respectively; the middle solid line expresses the phase boundary of 
spin component with Zeeman level m = 0, and keeps invariant when p = 0.08zt; the five solid lines 
represent the phase boundaries of spin components with Zeeman levels m = —2, —1, 0, 1, 2 respectively 
(from interior to external); the four dashed lines represent the phase boundaries of spin components 
with Zeeman levels m = —2, —1, 1, 2 respectively (from interior to external). 
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